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ABSTRACT 

In a multiwavelength program dedicated to identifying optical counterparts of faint persistent 
X-ray sources in the Galactic Bulge, we find an accurate X-ray position of SAX J1712. 6-3739 
through Chandra observations, and discover its faint optical counterpart using our data from 
EFOSC2 on the ESO 3.6m telescope. We find this source to be a highly extincted neutron star 
LMXB with blue optical colours. We serendipitously discover a relatively bright and large 
bow shock shaped nebula in our deep narrowband H a imaging, most likely associated with 
the X-ray binary. A nebula like this has never been observed before in association with a 
LMXB, and as such provides a unique laboratory to study the energetics of accretion and 
jets. We put forward different models to explain the possible ways the LMXB may form this 
nebulosity, and outline how they can be confirmed observationally. 
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1 INTRODUCTION 

The bright - X-ray luminosities in the order of 10^'' erg s"' and 
higher - accreting black hole and neutron star X-ray binary systems 
are amongst the best studied Galactic X-ray sources. Many of these 
show a rich phenomenology, displaying spectral state changes, dra- 
matic accretion outbursts and complex temporal behaviour on a 
wide range of temporal and spectral scales. However, an interesting 
sub-class of X-ray sources exists that are persistent in their X-ray 
behaviour, but have low X-ray luminosities. The monitoring pro- 
gram of the Galactic bulge source population with the RXTE satel- 
lite through PC A scans started in 1999 and has revealed many new 
bright X-ray transients. It also detected a large number of these 
faint, persistent sources, for which relatively little is known. 

We have undertaken a small survey using Chandra and the 
ESO 3.6m telescope to accurately localize a sample of 20 of these 
RXTE faint persistent X-ray source sources, and determine their 
nature: this class of sources is thought to harbour a wide variety 
of exotic types of systems (e.g. the class of the very short binary 
period sources - ultracompact X-ray binaries). The main results of 
this survey will be presented in a forthcoming paper (Wiersema et 
al in prep; hereafter W09). 

One of the sources of this survey was SAX J1712. 6-3739, 
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a candidate ultra-compact X-ray binary (UCXB: lln 't Zand et al.l 
[2007). In the case of a UCXB the mass accretion rate or luminosity 
may be a good indicator for the mass transfer rate from the donor 
star, and such low accretion rates suggest the possibility of degener- 
ate donors such a white dwarfs, making UCXBs important sources 
to study from a binary evolution point of view. However, their (op- 
tical) faintness makes UCXBs parti cularly hard to find an d unam- 
biguously identify (see for example [Nelemans et al.ll2006h . 

At optical wavelenghts, several components contribute to the 
light received from a LMXB, e.g. the accretion disk, the secondary 
star and a possible jet component. To broadly separate these vari- 
ous contributions and eliminate other, perhaps less likely, possible 
source identifications (e.g. accreting white dwarfs, high mass X-ray 
binaries; see W09 for statistics from our survey), we observed the 
sources in the survey of W09 in four optical bands. 

In our narrowband H tt imaging of SAX 11712.6-3739, we dis- 
cover a large, bright, roughly bow shock-shaped nebulosity, with 
an origin consistent with the position of the Chandra counterpart 
of SAX 11712.6-3739, and unlike any other previously discov- 
ered nebulosity associated with an X-ray binary. The importance 
of this discovery is clear: the production of the H a nebula in SAX 
11712.6-3739 is undoubtedly related to the accretion produced pho- 
ton flux or the mechanical work of an accretion-associated jet slam- 
ming into the ISM (and the kinetic energy of the binary), and in 
both cases the morphology and flux of the nebulosity can be used 
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Figure 1. The PCA bulge scan lightcurve spanning several years, taken 
from the PCA Bulge scan web database. Indicated by dashed lines are the 
times of our Chandra observation and the two XMM-Newton slews. 



to constrain the current and past energy output of this candidate 
neutron star UCXB. 



2 X-RAY OBSERVATIONS 

As part of the survey described in W09 we observed the BeppoSAX 
position of SAX J1712.6-3739 with the High Resolution Camera 
(HRC-I) onboard Chandra, on 2007 October 23, with an effective 
on-source exposure time of 1.176 ks. The data were processed us- 
ing the CIAO tools (version 4.0) and standard Chandra analysis 
threads. Using the wavdetect routine, a bright point source is found 
with count rate 5.16 counts s"', at position RA: 258.15319 degrees, 
Dec: -31 .MAlTi degrees (J2000), with an astrometric uncertainty 
of 0.6 arcseconds (90% containment, including systematics). This 
position is within the errorcircle of the ROSAT position of SAX 
J17 12.6-3739, the best previous position of this source. Based on 
position and observed count rate we conclude that the Chandra 
source can be identified with SAX J1712. 6-3739. 

W e find that the XMM-Newton slew survey jSaxton et al.l 
l2008h lists a bright source consistent with the Chandra position: 
XMMSLl J171236.3-373841. The source is detected in two slews: 
the first during 2006 March 8 and the second during 2008 February 
27, at position RA: 258.15194 degrees, Dec: -37.644827 degrees 
(J2000; ~ 6" error). XMM-Newton EPIC-pn 0.2-12 keV count rates 
were high: 23 counts s"' in slew 1, and 64 counts s"' in slew 2, with 
elTective exposure times of 3.7 and 4.2 seconds, respectively. PCA 
observations at similar times also show a count rate difference of 
a factor ~ 3.5 (Figure [TJ, strengthening the identification of the 
XM M-Newton source wi th the PCA source even more. 

ICocchi et all ( 1200 ih use the BeppoSAX data to find N^ ~ \ - 
2 X 10^^ cmr^, and r ~ 2.2. We extract spectra from the XMM- 
Newton slew data, and fit them with an absorbed powerlaw, where 
we fix r = 2.2 in order to get a good estimate of the absorption, 
needed to interpret the optical data (see Section|4l(, as there are not 
sufficient counts to fit for both parameters. We find Wh = \.i9*Q^^^x 
10^^ cm-2 (90%) for slew 1, and Nn = 2.36!;!;5^ x 10^^ ^m^^ (90%) 
for slew 2. This is consistent w ith, though on the high end of, the 
findings of lCocchi et all j200lh . 



3 OPTICAL OBSERVATIONS 

Observations were performed in Service Mode with the EFOSC2 
instrument on the ESQ 3.6m telescope at La Silla, Chile, at the 
BeppoSAX position of SAX J1712. 6-3739, as part of the survey 
described in W09. This configuration gives a field of view of 
~ 4.5 X 4.5', with a pixel scale of 0.31" per pixel on sky (2x2 
binning). Data were taken on 2007 July 1, in the following filters: 




Figure 2. EFOSC2 H a image of the field of SAX J1712.6-3739, excluding 
the areas on the sides that are most affected by vignetting. The position of 
the bright Chandra source is indicated with a white circle. Note that this cii- 
cle is much larger than the Chandra errorcircle, for presentation purposes. 
The inset gives a rough guide to the main features. 



Bessel V (ESO filter #641; total exposure time 10 minutes), Bessel 
R (#642; 10 minutes), Gunn i (#705; 10 minutes) and narrowband 
Ha (#692, 15 minutes). Conditions cannot be considered photo- 
metric with certainty, and as such an accurate absolute calibration 
using photometric standard stars cannot be performed. A full log of 
the optical observations can be found in Table[T] 

All data were overscan corrected, trimmed, flatfielded and cos- 
mic ray cleaned with the usual procedures, using tasks in IRAF. 
Data were taken close to the Moon, giving rather high sky back- 
ground, resulting in somewhat lower limiting magnitudes for these 
bands than can be expected under the seeing and airmass conditions 
given in Table [T] 

WCS calibration of the optical frames is carried out by match- 
ing sources in the frames to USNO-B2 source lists. The resulting 
WCS solutions (RMS 0.3") were subsequently refined using the 
Two Micron All Sky Survey (2MASS) sources detected in the field. 
Resulting RMS values of the best astrometric fit for each filter are 
0.13" for the R and i band frames (using 590 2MASS stars), and 
0.16" for the V and Ha frames (400 stars). 

As this field is in the Galactic plane (latitude b = 0.92 degrees) 
the field is crowded, particularly in i and R band. We therefore used 
IRAF DAOPHOT tasks to perform iterative point spread function 
(PSF) fitting photometry: after detecting a bright sample of stars 
with DAOFIND, we fit a first order spatially variable empirical PSF 
to bright stars (using several iterations to subtract off neighbour 
stars, refine PSF star selection and use progressively larger fit radii). 
After PSF photometry using ALLSTAR we subtract this list of stars 
and iterate with a lower detection threshold. 

EFOSC2 is a focal reducer instrument, and as such suffers 
from a small amount of vignetting, and the PSF for sources near 
the chip edges is distorted. The position of the optical counterpart 
of SAX J1712. 6-3739 is sufficiently away from the chip edges that 
this does not affect our results (see Figure |2}. 
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Table 1. Log of our observations of SAX J1712.6-3739. The 5ct limiting 
magnitude listed is derived for the entire EFOSC2 frame, and is therefore 
not equivalent to a detection limit of a specific source on a specific place on 
the chip (i.e. when confusion / blending plays a role). Seeing is measured 
from the PSF on the frames, excluding the areas affected by vignetting. 





Band 


Observing time 
start (UT) 


Exp. time 
(s) 


Seeing 
(arcsec) 


Airmass 


Lim. mag. 

(5 IT) 


EF0SC2 


R 


01/07/2007 04:28 


2x300 


0.65 


1.04 


23.4 




V 


01/07/2007 04:39 


2x300 


0.76 


1.06 


23.5 




I 


01/07/2007 05:51 


2x300 


0.85 


1.07 


23.0 




Hff 


01/07/2007 05:02 


900 


0.71 


1.08 




Chandra 


HRC-I 


23/10/2007 


1.176 ks 









Table 2. Results on the optical counterparts of SAX J1712.6-3739. The op- 
tical positions are measured on the I band frame, where S/N is greatest, and 
include the absolute WCS calibration uncertainty to the 2MASS astrometry 
(see Section[3). Please note that the night was likely not photometric, and 
therefore the errors given on the magnitudes below do not include an eiTor 
from absolute photometiic calibration. 





Position (J2000) 


Error (") 


I 


R 


V 


SI 
S2 


258.15318-37.64501 
258.15376-37.64465 


0.25 
0.25 


21.45 + 0.06 
21.38 ± 0.04 


22.61 +0.12 
23.05 ±0.13 


23.93 + 0.26 
> 24.1 (3a-) 



To derive approximate magnitudes, we use EF0SC2 zero- 
points and La Silla average atmospheric extinction terms from the 
EFOSC2 Quality Control webpages, calibrating the Gunn i band 
data to I. All magnitudes are given in Vega magnitudes. While the 
absolute photometry is likely not highly accurate, the differential 
colours of objects in the field should be reasonably reliable, allow- 
ing a classification of sources. Note that due to the low Galactic 
latitude no SDSS photometry is available for field star calibration. 

We use the V, R, I and Ha photometry lists to construct 
colour-magnitude diagrams (CMD) and colour-colour diagrams 
(CCD). For each colour magnitude diagram point sources are 
matched for the relevant filters only. The R - Ha colour cannot 
be reliably quantified, as no spectrophotometric standard was ob- 
served. However, most sources form a narrow strip of R - Ho' 
colour (i.e. no unusally strong H a emission or absorption with re- 
spect to the R band continuum), which we shift to an average value 
of R - H a ~ 0, to guide the eye. As the H a line is not situated 
directly in the center of the Bessel R band, we use a R - i vs R - 
Ha diagram whenever possible, following the procedure outlined 
in Witham et al. (2006) to identify Ho- emitters near the Chandra 
position. To find potential optical counterparts to the X-ray source, 
we add the (90%) RMS of the WCS calibration of the optical im- 
ages to 2MASS in quadrature to the (90% containment) X-ray posi- 
tion error for each filter, and search for optical point sources within 
the resulting error circle (Figure 3). 



4 SAX J1712.6-3739 

This source has been discovered in 1999 (in 't Zand et al 1999), 
and is a persistent source, detected by ROSAT, BeppoSAX and 
INTEGRAL (alternative identifiers INTEGRALX 47 and IRXS 
J171237. 1-373834). It has been classified as a neutron star LMXB, 
on the basis of several type I X -ray bursts, which set the distance 
at 6 - 8 kpc iCocchietaUboOll ; see e.g. in 't Zand et al. 2007 for a 
discussion on past observations). 



3 




Figure 3. A small section of the EF0SC2 i band image of the field of SAX 
J1712. 6-3739. The position of the bright Chandra X-ray source is indicated 
with a white circle. Two optical sources are located close to the Chandra po- 
sition, labelled SI and S2. SI is located significantly closer to the Chandra 
position and has bluer colours than S2. 




R - 1 V-R 

Figure 4. CMDs of the field of SAX J17I2.6-3739, with the optical candi- 
dates SI indicated in red, and S2 in green. Errorbars are omitted for clarity. 
S 1 is clearly a blue source, and is one of the faintest R band sources with a 
V band detection. S2 is not detected in V. 

The EF0SC2 H a data shows patchy nebulosity, which is con- 
firmed through much larger field of view Ha d ata retrieved from 
the SuperCosmos repository jParker et alhoOSi) : the X-ray source 
is located in an area abundant in bright difi^use clouds of gas. In the 
EF0SC2 data a bright nebular source is obvious, which appears 
distinctly jet or bow shock shaped and spans most of the H a im- 
age (see Figure |2j. The R band data are, not surprisingly, the only 
wideband data showing the central part of this nebula, consistent 
with an (H a) line emission dominated nebula. 

The optical data shows two point sources very close to the 



4 K. Wiersema et al. 



Chandra position, both detected in R and i, see Figure[3]and Table 
|2] The more southern source (hereafter called SI) is marginally de- 
tected in V (~ 4(t), somewhat hindered by the proximity to diffrac- 
tion spikes of a brighter star. The other source (hereafter called S2), 
to the East of the Chandra position is not detected in V. The posi- 
tion of SI (see Table |2j is close to the X-ray position (1.0" from 
the center of the Chandra errorcircle, i.e. consistent within 2a), 
whereas S2 is considerably further away (2.1"). 

At the positions of these two sou rces no point source i s found 
in the H a frame. Using the method bv lWitham et alj j2006l) to plot 
sources inaR-HavsR-I diagram to select H a emitters, we find 
that the H a detection limits are not constraining for either source. 
SI has clearly very blue (intrinsic) colours, see Figure]?] consistent 
with expected colours of a LMXB (though note that this depends 
on the adopted extinction value, which is high and very uncertain). 
Position and colours make SI a more likely counterpart - though 
spectroscopy or photometric variability will be required to confirm 
the association. 

lln 't Zand et alj ( l2007h suggest from a literature data compila- 
tion that this source is a candidate UCX B. They derive a reddening 
by using the IPredehl & Schmit^ i 19951) rel ation betw een Galactic 
A^H and Ay, using A^h = 1.3 x 10^^ cm"^ jCocchi e t al. 2001 and 
lln 't Zand et al.l2007l a value roughly consistent with our estimate, 
see Section |2]l and find Ay = 7.3. In combination with the dis- 
tance of ~ 7 kpc (Section |4). lln 't Zand et all j2007l) predict a vi- 
sual magnitude of ~ 26 for a UCXB and ~ 22 if the source is a 
non-ultracompact binary. Sources SI and S2 have V magnitudes 
(or limit in the case of S2) in between these two values, see Table 
|2l but obviously this is rather dependent on the uncertain extinction 
towards this source. Using this same extinction estimate and dis- 
tance we find My — h2.4 for SI. While this is faint, it is somewhat 
brighter than the absolute magnitudes of UCXB candidates, and 
comparable to standard short period LMXB values (van Paradijs & 
McClintock 1994), though we note the uncertain optical extinction. 
Time resolved optical spectroscopy of this system would be able to 
reveal the binary period and therefore the UCXB nature. 



5 THE NATURE OF THE NEBULA 

The Chandra position coincides with the position from which the 
jet- or bow shock shaped nebula visible in H a appears to originate, 
see Figure (2] strengthening a connection between this object and 
the nebulosity. Specifically, an extension of the two approximately 
straight band / stripes of H a emission meet near the location of the 
LMXB. The probability of chance coincidence of a LMXB with 
an, unassociated, bow shock shaped nebula can however not be ro- 
bustly calculated, as it would require searching H a surveys for sim- 
ilar features unrelated to LMXBs, while the definition of "similar" 
is necessarily influenced by the detection of this current, unique, 
source. It is important to emphasize here that a theoretical mode l 
predicting a system similar to this observed one jHeinz et al.l2008l) 
has been made prior to our discovery of this system. 

We may expect the LMXB to move roughly with the Galac- 
tic plane, as a large kick is required to make it move away from 
the plane. The projected direction of motion of the LMXB inferred 
from the bow shock shape is towards the Galactic plane, i.e. moving 
towards lower Galactic latitude, consistent with this picture. 

In the following we continue on the premise that the nebula 
is indeed produced by the LMXB, and discuss the two most likely 
mechanisms to produce this nebula, discussing also which future 
observations may help distinguish these models. 



5.1 Morphology 

The nebula seems to exist of two nearly parallel trails of H a emis- 
sion that appear to originate close to the LMXB, and a diffuse, near- 
semicircular region of emission just 'in front' (to the south-east) of 
SAX 11712.6-3739. The distance from the LMXB to the leading 
edge of this diffuse region is ~ 30", which corresponds to a distance 
of ~ 3 X 10'* cm (adopting d = 1 kpc). The trails appear to fade 
and widen with distance from the LMXB except initially; the max- 
imum surface brightness of the trails is located at a distance ~ 20" 
and ~ 40" from the LMXB for the northern and southern trails, 
respectively. The width of each trail is ~ 9" close to the LMXB, 
widening to ~ 33" at a distance of 100", for the northern trail. 



5.2 A jet-powered nebula? 

Accretion is on almost all masses and timescales observed to re- 
late to the formation of high velocity jets, from the extremely en- 
ergetic gamma-ray bursts to the more sedate outflows of young 
stellar objects. Evidence for jets in X-ray binaries is found from 
direct detection of fast moving radio or X-ray emis sion from the 
jets themselves or shocks downstream (e.g. iTingav et al, 1 19951 : 
iFomalont et al.ll200ll : lcorbel et alfcOoH: iTudose et al II2008I) , from 



transient (unresolved) radio emission lFendedl200l'l) . from spec- 
tral/polarimetric/timing studies in the infrared and optical (e.g . 
Corbel & Fe nder 2002; Malzac et al. 20041 iRussell et all 12007 j; 
Shahbaz et al.,,2008, ; .Russell & Fender ,20081 ) or in one case from 



Dopp ler-shifted emission lines from baryonic jets jMargon et al.l 
Il979l) . Recently, bow shock nebulae powered by the jets of Cyg 
X-1 and LMC X-1 have been identified, and their pro perties have 
been used to infer the time-averaged power o f the jets dGallo et al.l 
l2005l ; lRussell et al.l2007bl ; ICooke et al.l2007 ). Since SAX J1712.6- 
3739 has a hard spectrum and a steady luminosity we may assume 
it to have compact jets, detectable at radio frequencies, with typical 
flux densities ~ 0.1-1 miy (Migliari & Fender 2006), unfortu- 
nately fainter than the limit of the NVSS survey locally. 

ludging by the density of Ho- clouds in the field (from the 
SuperCosmos data), this source is situated in a dense interstellar 
medium. The long term interaction of a jet from a LMXB with a 
dense medium coupled with a high space velocity may give rise to 
bow shock-shape d trails of collisionally excited gas (see for details 
iHeinz et al.l2008l) . The requirement of a high space velocity means 
that LMXBs should more easily form jet trails than HMXBs. How- 
ever, sources with a morphology like this have not been observed 
to date, probably because for m ost LMXBs the local ISM density is 
not high enough ( lHeinzll2002h . The Heinz et al. (2008) model pre- 
dicts a morphology and size remarkably similar to the observed H a 
nebula (FigureO which would make SAX 11712.6-3739 the first 
confirmed jet trail source. In this scenario, we interpret the nebula 
as an edge-brightened cone-shaped bow shock (a 'head-tail source', 
like those seen associated with AGN) powered by two continuously 
replenished jets. The space velocity of the LMXB would exceed 
the velocity of the bow shock, forming trails behind the LMXB as 
opposed to shell-l ike stationary cocoons like the one of Cyg X-1 
dGallo et al]|2005l) . The emission region 'in front' of the LMXB 
is expected in this scenario, ionized by the approaching radiation 
field. The H a emission originates in this case from recombination 
of the shock-compressed gas near the shock front. Radio plasma 
consisting of relativistic particles and magnetic fields should fill the 
region behind the bow shock front - these trails may be detectable, 
as might radio-emitting bremsstrahlung radiation from the shell. 
The detection of such structures will confirm the jet-powered na- 
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ture of the nebula. When interpreting the nebula as a jet blown trail, 
the width of the trail provides importa nt constraint s on th e system. 
In particular, using equation 12 from Heinz et al. ( 2008) we find 
that the average jet power W is likely well below 10^' erg s"', not 
unexpected for a neutron star LMXB (though note that prototypical 
neutron star jet source Cir X -1 has a jet power close to 10^^' erg s"' ; 
see e.g. iTudose et al.ll2008l and references therein). In particular, 
the persistently active, neutron star UCXB 4U 0614+09 has been 
shown to produce powerful jets iMigliari et al 12006.) with jet power 
~ 10 ^^ erg s~', dominating the radio to nIR emission from the sys- 
tem iMigliari et al.ll2006l ; [vrtileldl2008 l, Migliari et al in prep). The 
jet powered scenario is therefore not unexpected. 

Ho- bow shocks with broadly similar morphology are occa- 
sionally found around fast moving isolated neutron stars, though 
much fainter and smaller than this source (e.g. lChatteriee & Corded 
l2002l ). powered by pulsar winds. Accretion inhibits strong pulsar 
winds in LMXBs but instead, accretion disc winds powered by ir- 
radiation on the disc could be powerful. However, disc winds are 
likely to be orders of magnitude weaker than the jets (see the case 
of 4U 0614-1-09 above), making it much more likely that the trails 
are powered by jets than by a disc wind. 



5.3 An X-ray-ionized nebula? 

Since SAX JI7I2. 6-3739 is a persistent X-ray source, local neb- 
ulosity will be photo-ionized by its X-ray and UV radiation field. 
The nebula surrounding LMC X-1 is partly ioniz ed in this fash- 
ion jPakull & AngebaultlfTgsi : ICooke et alll2007h . The morphol- 
ogy of the SAX J17 12.6-3739 nebula may be consistent with this 
interpretation if the LMXB has a high space velocity. However, the 
observed edge brightened morphology is hard to achieve in this sce- 
nario, while expected in the jet trail model. Furthermore we would 
need either a definite density gradient along the binary trajectory or 
a declining ionizing photon flux as a function of time to create the 
observed nebula morphology. 

Emission line ratios are commonly used as diagnostics of 
shock-ionized or photo-ionized gas (e.g. ICox & Ravmondll98^ : 
lAllen et alj|2008l) : for example the [S II] / H a line ratio is an in- 
dicator of shocked gas and the He II / H ratio of photo-ionized 
gas. Follow-up spectroscopy (or narrow-band imaging) of the trails 
will therefore constrain the nature of the nebula and probe the ve- 
locity, density and composition of the gas, constraining either the 
time-averaged power of the jets or the soft X-ray/EUV radiation 
field of the LMXB. Note that in both models discussed above the re- 
quirement of a high space velocity suggests a relatively high proper 
motion, which may be detectable with future instrumentation. 



6 CONCLUSIONS 

In a survey of faint persistent X-ray sources with Chandra and 
EF0SC2 on the ESQ 3.6m telescope, we find an accurate X-ray 
position of the neutron star low mass X-ray binary SAX J1712.6- 
3739 and identify its most likely optical counterpart: a faint source 
with blue colours with respect to field stars. Using the best known 
extinction and distance estimates, we find that a UCXB nature of 
SAX JI7I2. 6-3739 can be accomodated by the measured magni- 
tudes. In the H a data we serendipitously detect a large nebula with 
a distinct bow shock shape, likely associated with the LMXB. The 
nebula is likely formed through either the interaction of a energetic 
jet from the LMXB with the interstellar medium, or through in- 
tense UV/X-ray radiation. In both cases a high space velocity is 



required, and in both cases the nebula can provide unique insight 
into the energetics of this system. 
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